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1. Introduction 

The Nambu-Jona-Lasinio (NJL) model is an 
effective theory of QCD at intermediate energies 
where pointlike interactions between quarks re- 
place the full gluon mediated description of quark 
interactions. The model has been widely used 
to study cold dense quark matter (QM) 0, be- 
cause in this region of the QCD phase diagram 
the explicit gluonic effects are expected to be mi- 
nor. Furthermore, recent developments of the 
NJL model have shown the possibility of a re- 
alistic description of single nucleons and stable 
nuclear matter (NM) 0. Hadronization tech- 
niques 0] in principle enable us to study both 
NM and QM within the framework of a single 
model, which also accounts for the internal quark 
structure of the free nucleon [5]. The purpose of 
this letter is to investigate the phase diagram of 
this model, extending the work of ref. 0] to the 
case of isospin asymmetric matter, and to present 
results for the equation of state (EOS) of charge 
neutral matter and the structure of compact stars 
□ 
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Many recent theoretical studies have suggested 
that matter at high densities and low tempera- 
tures is in the color superconducting quark matter 
(SQM) phase |7J. The number of flavors present 
will depend on the effective quark masses in QM. 
In particular the behavior of the strange quark 
mass at high density, which is not well known, 
plays a critical role in determining the structure 
of the favored QM phase In NJL- type models, 
the strange quark turns out to be heavy enough 
to favor a transition to the 2-flavor color SQM 
state 0, whereas the 3- flavor state occurs at still 
higher densities. There are many recent investi- 
gations on the dependence on the strange quark 
mass [911011 1] which tend to support this sce- 
nario, provided that the quark pairing interaction 
is strong enough. 

In the present work we also examine the ef- 
fects of color superconductivity on the phase di- 
agrams and the EOS at high densities. It is a 
first attempt toward the goal of describing pos- 
sible mixed NM/SQM phases in one framework. 
We will use a description which, in normal NM, 
avoids unphysical thresholds for the decay of the 
nucleon into quarks, thereby simulating the ef- 
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feet of confinement 3 . We will show the result- 
ing EOS including mixed phases, and discuss the 
consequences for compact stars. 

We derive both NM and QM phases from the 
flavor SU(2) NJL-Lagrangian, 

L = ^ - m)i> + Gjr ((^) 2 - (^ 75 tV) 2 ) 

+ G s (^ 7 5Cr 2 /3' 4 ^ T )(^ T C- 1 75T2^V) , (1.1) 

where we show only the interaction terms rele- 
vant for our discussions 1 . Here to is the current 
quark mass, ij) is the flavor SU(2) quark field, and 
the coupling constants G n , G w and G p character- 
ize the qq interactions in the scalar, pseudoscalar 
and vector meson channels, while G s refers to the 
interaction in the scalar diquark channel. 

The model is further specified by a regulariza- 
tion scheme, for which we use the proper-time 
scheme ^31 m this work. It is characterized by 
an infrared cut-off (Air) in addition to the usual 
ultraviolet one (Auv)- In the vacuum, the pa- 
rameters of the model are determined as follows: 
We fix Air = 200 MeV, and choose Auv, m and 
Gjr so as to reproduce /„• = 93 MeV, m, = 140 
MeV, and constituent quark mass Mq = 400 MeV 
via the gap equation at zero density. We will set 
M = in QM because the quark mass is already 
very small in the region where the transition to 
QM occurs. We have found that including the ef- 
fective quark mass in this phase does not change 
the structure of the phase diagrams and has little 
effect on the EOS. 

2. Nuclear matter 

The nucleon is constructed as a quark-diquark 
bound state jS], making use of the scalar diquark 
interaction term in and the Bethe-Salpeter 
equation to get the scalar diquark mass, M s . 
The interaction with the spectator quark is de- 
scribed by the quark exchange (Faddeev) kernel, 
for which we use a momentum-independent ap- 



1 We note that every 4-fermi interaction Lagrangian can be 
decomposed, as in Eg. II. II . into various qq and qq chan- 
nels by using Fierz transformations 5;. In the last term 
of Eq.JTTJ, C = J7270 and f3 A = y r Zj2\ A (A = 2,5,7) 
are the color 3 matrices. 



proximation in the finite density calculations re- 
ported in this letter. (This corresponds to the 
"static approximation" of the Faddeev kernel - 
see Refs. |15l3] for details.) The coupling constant 
G s is chosen to reproduce the free nucleon mass, 
M N0 = 940 MeV. 

In the mean field approximation, the NM phase 
is characterized by composite neutrons and pro- 
tons, moving in scalar and vector mean fields. 
There is also a non-interacting sea of electrons 
in chemical equilibrium with the nucleons (fi e = 
— [J> p , where \i denotes the chemical poten- 
tial) 2 . The form of the effective potential in 
the mean field approximation has been derived 
for symmetric NM in Ref.flj starting from the 
quark Lagrangian and using the hadroniza- 
tion method. This can be easily extended to the 
isospin asymmetric case. It has the form 



y(NM) = + y N 



4G. 



where 
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V N = -2 £ / 



d 3 k 
J2nf 



a— p,n 

90£ (/£ ~E N (k)) (2.3) 

describes the Fermi motion of nucleons mov- 
ing in the scalar and vector mean fields. We 
used E N (k) = ^ M 2 , + k 2 , where M N {M) is 
the nucleon mass in-medium, which is obtained 
from the pole of the quark-diquark T-matrix. 
The effective chemical potentials are defined as 
/i* = fi a — 3wo T 3/9o in terms of the mean vec- 
tor fields to = 2G W (NM|^V|NM) and p° = 
2G p (NM|V' t T3^|NM). 

2 At some density the components may change, as neg- 
atively charged kaon or pion condensates |lti| as well as 
muons replace the electrons and more massive hadron 
species replace nucleons 1171 . but here we confine ourselves 
to the most simple picture. 
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The constituent quark mass, M, and the mean 
vector fields in NM are determined by minimizing 
y(NM) £ or £ xe( j chemical potentials fi p and fi n . 
The parameter G u is fixed by the requirement 
that the binding energy per nucleon of symmetric 
NM passes through the empirical saturation point 
(baryon density pa = 0.16 fm -3 and Eb/A = 15 
MeV) 3 , and the parameter G p is adjusted to the 
empirical symmetry energy (04 = 32 MeV at po). 
The resulting values of the parameters are shown 
in the column "NM" of Table 1. 

As compared to chiral models for point- 
nucleons, the important property which leads 
to saturation of the NM binding energy in this 
approach is the positive curvature of the func- 
tion Mjv(M), which reflects the internal quark 
structure of the nucleon and which works effi- 
ciently only if there are no unphysical thresholds 
for the decay of the nucleon into quarks 3 . In 
our method, these thresholds are avoided by the 
choice of the proper-time regularization scheme 
with A m > [H]. 





NM 


QM 


m [MeV] 


16.93 


17.08 


G^GcV- 2 ] 


19.60 


19.76 


AuvpvleV] 


638.5 


636.7 


Am [MeV] 


200.0 





r uj = G u 1 G v 


0.37 





r p = Gp/G-x 


0.092 





r s = Gs/Gk 


0.51 


free parameter 



Table 1 

Parameters used for nuclear matter (left column) 
and for quark matter (right column). The proper 
time regularization scheme is used in both cases. 
Because Air is set to zero in the QM case, the 
parameters m, G^ and Atjv differ slightly from 
the NM values in order to obtain f„ = 93 MeV, 
= 140 MeV, and M Q = 400 MeV. 



3 We recall from Ref.|3 that is the only free param- 
eter for NM. With the limited number of parameters in 
this simple model it is not possible to ensure that the cal- 
culated binding energy curve also has a minimum at the 
empirical saturation point. Instead it occurs at p = 0.22 
fm" 3 and E B /A = 17 MeV. 



3. Quark matter 

The effective potential for QM in the mean field 
approximation, allowing for the possibility of di- 
quark condensation, has the form 

V« M > = Kac + V Q + V A -^5 (3.1) 

where the vacuum part, V v&Cl is given by (|2.2|1 
and 

a— u,d 

e(fi a - E Q (k))(n a -E Q (k)) (3.2) 

describes the Fermi motion of quarks with 
chemical potentials fi u and pa, and Eg(k) = 
V M 2 + k 2 . The term Va describes the effect of 
the pairing gap and is given by 

+ ln k 2 -(E a -^) 2 \ + 6G S ' W 

where e±(k) = (E(k) ± /J q ) 2 + A 2 and E± = 
\E(k) ± (i q \. Here we used the isoscalar and 
isovector combinations fi q — (/i u + /id)/2 and 

M/ = (M« - Md)/2 4 . 

A detailed discussion of vector meson poles in 
Ref.0] has shown that the vector- type interac- 
tions should be set to zero in QM. In addition, 
we set Air = 0, as the infrared cut-off simulates 
confinement effects which are not appropriate in 
QM. Also, as stated in the Introduction, we will 
set M = in QM, since the effects of quark mass 
are small when constructing the phase diagrams 
and the charge neutral EOS. The resulting pa- 
rameters are shown in the column "QM" of Ta- 
ble 1 . In the main part of this paper we treat G s 
as a free parameter in QM to investigate the de- 
pendence on the pairing strength, and comment 
on the important question of consistency with the 

4 We mention that in principle one needs a further chemical 
potential for color neutrality (fig) in QM. However, for 
the 2-flavor case pig turns out to be very small 2 23 , and 
therefore we neglect it here for simplicity. 
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value derived from the free nucleon mass at the 
end of Sect. 4. 

The gap A is determined by minimizing the 
effective potential for fixed quark chemical po- 
tentials. In the following discussions, we will dis- 
tinguish the normal quark matter (NQM) phase, 
which is characterized by A = 0, from the SQM 
phase (A > 0). 

4. Results 

In order to construct the phase diagram, we 
compare the effective potentials in the NM, NQM 
and SQM phases for several fixed chemical poten- 
tials for baryon number and isospin (ps and /Ltj), 
which are related to the chemical potentials of the 
particle species by 

p a = p B ± M/ ( a = P, n ) , 
Ma = — ± M/ \ a = u > d ) > 

fi e = -2fii. (4.1) 

Fig.l shows the phase diagrams for several 
choices of the pairing strength in the SQM phase. 
The black, dark and light regions indicate the 
phases with the lowest effective potential, and 
the ± indicate the sign of the total charge den- 
sity 5 . The corresponding plots in the plane of 
baryon and charge density are shown in Fig. 2. 
The mixed phases appear as white regions in this 
Figure. The phase boundaries, which are single 
lines in Fig.l, appear as two lines facing each 
other in Fig. 2. By connecting the correspond- 
ing end points on the boundaries by straight lines 
(dashed lines in Fig. 2), we can divide the white 
region into sections which correspond to mixtures 
of the two phases facing each other. The p,j = 
axes in Fig.l correspond to the upper-most line 
running from NM to QM in Fig. 2, and because 
we consider only the case fxi < the upper left 
parts of the diagrams in Fig. 2 are left empty 6 . 

5 The density of particle A = p,n,u,d,e is obtained as 
PA = —dV/dfiA- In particular, the baryon and charge 
densities in NM are pg = p p + p n , pc = Pp — pe, and in 
QM p B = (p u + p d )/3, p c = 2/3p u - l/3p d - p e - 

6 In principle we could extend the phase diagrams to the 
region where pj < by admixing positrons instead of 



For each case in Fig.l, the charge neutral EOS 
corresponds to the line separating the positively 
and negatively charged regions. For example in 
the first diagram (r s = 0), the charge neutral 
EOS begins in the pure NM phase, then there 
is a mixed (+NM/-NQM) phase and finally a 
pure NQM phase. Each point on the bound- 
ary between the NM and NQM phases satis- 
fies the Gibbs conditions, since P^ nm \pb, Pi) = 
p(JVQM)( MBjM7 ). The EOS in the mixed neu- 
tral phase is found by using the method of Glen- 
denning PB], that is, the volume fraction of the 
NM phase is determined by the requirement of 
charge neutrality as x^ NM ^ = Pc M ^/(p c ® — 
pi NM) ), i« M ' = 1 - x^ NM \ Physically this 
method implies that the mixed phase begins with 
charge neutral NM, and then as the charge of 
NM becomes increasingly positive, regions of neg- 
atively charged NQM form, such that the mixed 
phase remains globally charge neutral 7 . The same 
sequence of neutral phases can also be seen in the 
first diagram of Fig. 2 by following the horizontal 
line pc = 0. 

As we increase the pairing strength in the QM 
phase, the regions where SQM is the ground state 
extend. When r s =0.1 we have a mixed (+NM/- 
NQM) phase, and then come to a triple point 
where all three phases meet. This is the same sit- 
uation as investigated in Ref. [23 , and leads to 
a region of constant pressure in the EOS, where 
all three phases are mixed. This is illustrated in 
Fig. 2 by the triangular region; that is, when the 
p c = line passes through this region all three 
phases are present. The volume fraction of SQM 
begins at zero on the left hand side of the trian- 
gle and increases while the volume fraction of NM 
decreases until it reaches zero on the right hand 
side of the triangle. The NQM phase occupies the 
remaining volume fraction, which varies continu- 
ously between the boundaries of the triangle. In 
this region the baryon density is increasing, while 
the pressure remains constant. However, within 
compact stars the pressure must always be de- 
electrons, but clearly the matter in this part of the phase 
diagram is always positively charged and is not relevant 
to the charge neutral EOS. 

7 For such a mixture one may calculate what sizes and 
shapes are favored for each component 1191 . 
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Figure 1. Phase diagrams in the plane of chemical potentials fis, fJ-i for baryon number and isospin for 
various choices of r s . The black, dark and light regions correspond to the NM, NQM and SQM phases, 
and the ± indicate the sign of the total charge density including the electrons in chemical equilibrium. 
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Figure 2. Phase diagrams in the plane of densities ps, pc for baryon number and charge for various 
choices of r s . The black, dark and light regions correspond to the NM, NQM and SQM phases, the white 
regions separated by the dashed lines correspond to the mixture of two phases, and the triangular region 
for the cases r s = 0.1 and 0.15 involves a mixture of three phases. The upper left regions in each diagram 
are left empty, because we consider only the case /xj < 0. The dotted line indicates charge neutrality. 




Figure 3. Left figure: Charge neutral EOS (pressure against baryon density) for pure NM (solid line), 
the transition to NQM (case r s = 0, dashed line), and the transition to SQM (case r s = 0.25, dotted 
line). Right figure: The corresponding star masses as functions of central density. 



creasing as a function of radial position. Thus 
the three component mixed phase cannot occupy 
any finite volume within a star and there will be 
a discontinuity in the star's density profile. 

If we further increase the pairing strength to 
r s = 0.15, a charge neutral SQM state be- 
comes possible. This transition from NM to 
SQM involves three intermediate mixed phases 
(NM -> +NM/-NQM -> +NM/-NQM/+SQM -> 
-NQM/+SQM ->■ SQM). Again, the three com- 
ponent mixed phase occurs at the triple point 
in Fig.l, which corresponds to the triangle in 
Fig. 2. At still larger pairing strengths the NQM 
phase becomes unfavorable and the EOS involves 
just NM and SQM. When r s = 0.25 we start 
with a neutral NM phase, and then enter a (- 
NM/+SQM) mixed phase before arriving at the 
neutral SQM phase. In Fig.l the line along the 
mixed phase in this case is only very short, which 
means that the pressure changes in the mixed 
phase are small. For (r s >0.3) the SQM phase 
almost completely expels the NM phase, which 
indicates an upper limit for the pairing strength 
which is consistent with the findings of Ref.0] for 



the isospin symmetric case. 

In Fig. 3 we show two examples of the charge 
neutral EOS, corresponding to the cases r s = 
and r s = 0.25, and the resulting star sequences 
obtained by integrating the TOV eauation|21|. 
(The results for the pure NM case are also shown 
for comparison.) In the first case, the system goes 
to the NM/NQM mixed phase at a baryon den- 
sity around 3.4po and the pressure of the mixed 
phase increases with increasing density. For cen- 
tral densities between 3.4po and 5.6po, stable hy- 
brid stars exist with a NM/NQM mixed phase in 
the center. For example, the maximum mass star 
has a radius of 11.6 km and the mixed phase is 
realized within r = 5.8 km. In the second case, 
the transition to the NM/SQM mixed phase oc- 
curs at around 2.3po ; the pressure in the mixed 
phase is almost constant, and around 4.6/?o the 
neutral SQM phase is reached. Stable hybrid 
stars exist with central densities between 2.3po 
and 3.4po with a very small region of NM/SQM 
mixed phase in the center and for central densi- 
ties between 5.6po and 13.4po stable quark stars 
may exist, which are composed of SQM in the 
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central region. For example, the maximum mass 
star for the case r s = 0.25 has a radius of 8.2 
km and the SQM phase is realized within r = 6.0 
km. These results are qualitatively similar to the 
ones reported in Ref. |9] . Comparison of two cases 
shown in Fig. 3 indicates that color superconduc- 
tivity certainly can have a major effect on the 
EOS. With increasing values of r s , both the den- 
sity and the pressure corresponding to the phase 
transition to QM are significantly reduced, the 
EOS becomes considerably softer, and the star 
mass decreases for given central density. 

It is interesting to note that for all cases which 
we studied and where we have a phase transition 
from neutral NM to neutral SQM (without ad- 
mixing the NQM phase), the pressure changes 
only very little in the NM/SQM mixed phase. 
Then the situation becomes similar to the one 
where a naive Maxwell construction is applied 
after imposing the neutrality condition for each 
phase separately. Including the quark mass in 
the QM phase does not change our conclusions. 
In the r s = 0.25 case (which has the lowest tran- 
sition density) the mass causes only a slight de- 
crease in the pressure (AP « lMeVfnT 3 ) in the 
region of the mixed phase. This is too small an 
effect to be seen on the scale of the phase dia- 
gram and certainly does not lead to any qualita- 
tive changes in our results. 

We note that most of the recent calculations 
on SQM have been performed by using the 3- 
momentum cut-off scheme and in general the gap 
is found to be around 100 MeV in the interme- 
diate density region [221 • Our work differs in 
that we have used the proper time regularization 
scheme, which leads to larger values of the gap 
(between 300 and 400 MeV in the relevant den- 
sity region). However, qualitatively the situation 
is similar to the cases of strong diquark coupling 
discussed in Ref. |9ll0| . 

Finally, we would like to come back to the im- 
portant question of consistency between the val- 
ues for G s used for the single nucleon and for the 
pairing strength in SQM (see Table 1). In this 
work, we fixed G s for the single nucleon by fitting 
the nucleon mass. However, we have to note that 
there are further attractive contributions to the 
nucleon mass, most importantly pion exchange 



24 25 26 and contributions of axial vector di- 
quarks [22], which would lead to smaller values 
of G s . Indeed, using the expressions given in 
Ref. 28 for the pion exchange contribution to the 
nucleon mass, we find that the nucleon mass can 
be reproduced with r s = 0.4(0.25), where the two 
numbers refer to the case without and with fur- 
ther inclusion of axial vector diquarks 8 . Thus the 
inclusion of pion exchange and axial vector di- 
quarks for the nucleon mass will allow a common 
value of G s for the single nucleon and SQM. The 
details of this calculation will be discussed else- 
where [2H] • 

5. Summary 

By applying a flavor SU(2) NJL model to both 
NM and QM phases, we have studied the phase 
diagram for isospin asymmetric matter at finite 
density. We emphasize that the model, and in 
particular the regularization scheme which we 
used, describes the single nucleon and the satu- 
ration of normal NM, and therefore forms a basis 
to investigate the EOS at higher densities. We 
found that, as we vary the pairing strength in 
QM, several scenarios are possible. The charge 
neutral EOS may make a transition to NQM or to 
SQM, via either one, two or three globally charge 
neutral mixed phases. These transitions begin at 
small enough densities (2.3 - 3.4po) that the QM 
phase, or at least the mixed phase, may occur 
inside neutron stars. 
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8 If axial vector diquarks are included r a is significantly 
decreased since there is alternative attraction through this 
channel. In this case G a (the axial coupling constant) is 
fixed using the Fadeev Equation for the A and setting 
M A = 1232 MeV. 
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